Gibbs attractor: a chaotic nearly Hamiltonian system, driven by external harmonic 

force 
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A chaotic autonomous Hamiltonian systems, perturbed by small damping and small external 
force, harmonically dependent on time, can acquire a strange attractor with properties similar to 
that of the canonical distribution - the Gibbs attractor. The evolution of the energy in such systems 
can be described as the energy diffusion. For the nonlinear PuUen - Edmonds oscillator with two 
degrees of freedom the properties of the Gibbs attractor and their dependence on parameters of the 
perturbation are studied both analytically and numerically. 

PACS numbers: 05.45.-a, 05.40.-a 



I. INTRODUCTION 

The Brownian motion of a particle of mass m in the 
static potential U{r) can be described by the system of 
Langevin equations 



mr + 7r + VU (r) = f (t) 



(1) 



where 7 is the viscous damping parameter and f (t) is the 
random force (white noise) with Gaussian distribution 
and mean values 

{f^ it)) = 0, (/. {t) , /, it')) = 2DS,,S {t - t') ■ (2) 

here Sij is the Kronecker symbol and 5{x) is the Dirac 
delta-function 1]. The damping term and the random 
force provide a phenomenological description of the in- 
teraction of the particle (in the potential) with its en- 
vironment, that is frequently called a heat bath. The 
dynamical system determined by Eqs. JQl with 7 = 
and D = will be called isolated. 

The ensemble of Brownian particles eventually will 
come to an equilibrium with the stationary probability 
distribution in the phase space given by the equation 



Q (p, r) Nexp 



H{p,ry 



e 



(3) 



where p = mr is the particle momentum, H (p, r) = 
p^/2m + U (r) is the Hamiltonian function of the iso- 
lated system, 9 = Df-f is the temperature, and N is 
the normalization constant. The distribution Eq. ^ is 
known as the canonical, or Gibbs, distribution; it serves 
as a central point of the equilibrium statistical physics 
0. From Eq. J^J follows the equilibrium energy distri- 
bution, 



Q {E) = N'^ [E) exp 



(4) 



where 



^{E) = J S{E-H (p, r)) dpdr 



(5) 



is the energy density of the phase volume on the energy 
surface H (p, r) = E. 

The canonical distribution Eq. (0) holds for the equi- 
librium state of the system Eq. ^ irrespectively of di- 
mensionality of the configuration space d and regardless 
of the nature of motion of the isolated system, be it reg- 
ular (periodic, quasiperiodic) or chaotic. 

Let's replace the random force f(t) by a regular one 
that depends on time harmonically; the equations of mo- 
tion will take the form 



mr + 71" + VJ7 (r) = F sin wt. 



(6) 
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If the motion of the isolated system is strongly chaotic, 
that is, nearly ergodic on the energy surfaces in a wide 
range of energy values, then we can expect that the mo- 
tion of weakly perturbed system (7 and F are small) will 
be chaotic too. In parallel with the physical picture of the 
Brownian motion given above, the external force Fsinwi 
will slowly change the energy of the system, whereas the 
damping 7 will provide a sink for the excessive energy, 
thus creating a possibility for the equilibrium. 

The dissipative non-autonomous system Eq. © may 
demonstrate the chaotic motion on a strange attractor, 
that has much in common with the canonical distribu- 
tion Eq. (j3J. For this reason it will be called the Gibbs 
attractor. 

The main purpose of this paper is to give an example 
of a system with the Gibbs attractor and to describe its 
main features. They mostly depend on the kinetics of 
energy exchange between the perturbation and the iso- 
lated system, that can be described as a process of energy 
diffusion. 

The problem that we have formulated is on a cross- 
road of several lines of research in nonlinear dynamics 
and non-stationary statistical physics. First, it is linked 
to the theory of chaos in non- autonomous Hamiltonian 
systems. In this theory the concept of energy (or action) 
diffusion is used for description of the infinite chaotic mo- 
tion above the threshold of its onset - in the models like 
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periodically kicked rotator Q or a hydrogen atom in the 
microwave field 0, Q . Unlike our problem, in these mod- 
els the strong periodic perturbation is the source that 
thrusts chaos on the system. Secondly, our problem is 
related to the theory of energy absorption by chaotic 
systems with parametric modulation, that is usually de- 
veloped (both in quantum and classical approaches) for 
models of billiards with varying form |3j 13 ■ As opposed 
to our problem, the main concern here is the study of slow 
variations, whereas we are interested in the perturbation 
with frequencies that are comparable to the characteristic 
frequencies of the unperturbed systems. Furthermore, in 
both aforementioned theories the damping is completely 
neglected. Thirdly, our problem has relation to the ac- 
tively developing theory of the Brownian motion under 
the influence of colored noise S, _9, IQ^ JJJ, since one 
may consider the harmonic perturbation as a limiting 
(monochromatic) case of the narrow-band noise. In con- 
trast to our problem, in this theory the unperturbed mo- 
tion is mainly one-dimensional, thus regular, and damp- 
ing is considered strong. Lastly, we note a connection 
to the issue of the effects of weak noise and damping on 
the Hamiltonian systems, that was discussed recently in 
the context of the problem of decay of metastable chaotic 
states [Hill. 

The rest of the text is organized as follows. In Sec. 
II the main equation for the energy diffusion is deduced 
in two ways. Sec. Ill consists of the description of the 
unperturbed model, the PuUcn - Edmonds nonlinear os- 
cillator, and main features of its chaotic motion. Sec. IV 
comprises the description of main properties of the Gibbs 
attractor in the model - its limits in the phase space, en- 
ergy distribution, conditions of existence and character 
of the energy correlation. Sec. V contains the concluding 
remarks. 



transitions is high, and the rate of these transitions can 
be described by the Fermi golden rule. With the account 
of dependence of matrix elements and density of states 
on the energy within the transition range, we can obtain 
for the rate of transitions the expression 



27r 
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S± 



fku 



S' + S^- 



(7) 



where S = Sx {E, lu) is the power spectrum of the co- 
ordinate X and p — p {E) is the density of states of the 
isolated system, both taken at the energy E, and dashes 
mean differentiating in energy p^ . 

The resonant absorption and emission of quanta pop- 
ulate narrow bands of levels, that are located on the en- 
ergy scale around values E' = E ± khio with integer k. 
We denote the probability of finding the system in such 
band of states as Q{E'). With the account of one-photon 
transitions we can write the balance equation 



dQ (E) 



dt 

-Q [E + hw) 



= -Q (E) [W+ + 

Q{E~hw) W. 



(8) 



Assuming Q{E,t) to be a smooth function of E, we 
can expand the arguments in the second and third terms 
in the RHS of Eq. © to the order of inclusive. This 
yields the equation in partial derivatives, that describes 
the energy diffusion in the Hamiltonian system under the 
infiuence of the external harmonic field: 



dQ d_ 
'dt^dE 



{AQ) 



dE \ dE 



0. 



(9) 



Here the coefficients of drift A and of diffusion D are 
given by the expressions 



II. ENERGY DIFFUSION 

We will study the evolution of the distribution of en- 
ergy values Q{E,t) for the Hamiltonian chaotic system, 
perturbed by small damping and small harmonic force. 
For the derivation of the equation of evolution for Q{E, t) 
we at first neglect the damping (we shall restore it later, 
in Subsec. C). 



A. The quantum approach 

The simplest approach is to start from the quantum 
model of the unperturbed system. Let's assume that at 
the initial moment the system is in a stationary state 
\n) with the energy E. The external harmonic force will 
induce transitions with absorption (-I-) and emission (— ) 
of the quanta fuv. If the motion of the classical system is 
chaotic and the power spectrum of the active coordinate 
X is continuous, then for its quantum counterpart in the 
quasiclassical case the density of final states with allowed 



A{E,u;) = -iu^F^S- 
2 p 



D{E,Lj) = -Lj^F'^S. (10) 



It is essential that in the quasiclassical case, when the 
density of states could be expressed by the Weyl formula 



P{E) 



<^(E) 



(11) 



the coefficients A and D do not depend on Planck's con- 
stant h. The analogous derivation of the equation for 
the energy diffusion in the conservative one-dimensional 
system under the inffuence of white noise was used in 



syst ( 
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B. The classical approach 



The equation Eq. Q does not depend on Planck's 
constant, but the condition of its applicability does. The 
Fermi golden rule is based on the assumption of the res- 
onant character of the transitions, which is justified only 
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when the transition rate W is much smaller then the 
perturbation frequency uj. Thus the diffusion coefhcient 
must obey the strong inequality 



D < h^oj^, 



(12) 



that is too restrictive for small H. 

However, we can rederive Eq. classically. In the 
zeroth approximation we neglect the influence of the ex- 
ternal force on the law of motion x{t). Then the in- 
stantaneous rate of the energy change is E = x{t)F{t). 
Variation of the energy (for the time interval T), A(r), 
in this approximation vanishes on the average for sym- 
metry reasons, (A(T)) ~ 0, whereas its averaged square 
is 

T T 

(A2(T))=y" j B,{h-t2)F{h)F{t2)dtidt2, (13) 



where B.^jir) is the correlation function of the x- 
component of velocity. If the unperturbed motion is er- 
godic, then the correlation function is determined by the 
microcanonical average 



1 



X (0) X {t)5 [E-H (p, r)) dp dr. 



(14) 

where x(0) and x{t) are taken on the trajectory that 
starts at t = at the phase point {p, r}. For times T 
much larger than the time of decay of velocity correla- 
tions Tc, we can rewrite Eq. H14|) in the form 



(A^ (T)) 



1 OO 

j dO j By (r) cos ujrdr. 



(15) 



The internal integral is proportional to the power 
spectrum of velocity, Sy{E,uj). Since Sy{E,uj) = 
uP'Sx {E, ui), we obtain for the coefficient of energy diffu- 
sion the expression D {E, lj) = {'k/2)uj^F'^Sx {E, lo), that 
coincides with the one given above in Eq. (|10|) . 

Let's assume that in the initial state the system has the 
energy Eq and denote Qo{E,t) the energy distribution 
with this initial condition, Qo{E,0) — 6{E — Eq). The 
equation Eq. H15|l can be rewritten as 



(A2(r)) =2 / D{Eo)de. 



(16) 



The distribution Qq spreads with time, including energy 
surfaces with different densities <& and different diffusion 
coefficients D. In the first approximation we can take 
this into account by averaging of these functions with 
the evolving distribution Qf){E,t): 



(A^ [t))=2 



<i>(i?o) 



dEQfi {E, 9) G {E) (17) 



where G (E) = ^ {E) D (E). In the time interval in which 
the distribution QQ{E,t) can be considered narrow, the 
function G{E) can be expanded in the Taylor series to 
the first order in A = £^ — £^o- Thus we come to the 
equation 



(A^ (t)) 



[D{Eo) 



G'(Eo) 
^{Eo) 



{^m]do. 



(18) 



If the system on the average changes its energy with a 
constant rate, (A {t)) = at, then (A^ (t)) = a^p _^ 2Dt. 
Substituting this expressions in Eq. H18|l . we obtain 



D = D{Eo), 



G' (Eq) 

^{Eo) ■ 



(19) 



The average rate of the energy variation in the state with 
given distribution Q{E, t) can be obtained from the equa- 
tion Eq. 10: 



{E) = J + D')Q {E, t) dE. 



(20) 



From Eqs. p9|l and ^2[)\ for the drift coefficient we obtain 
the expression A{E,u}) = D{E,uj) ($'(£')/$(£')), that 
coincides with the one given above in Eq. (|10|) . since 
the density of states is proportional to the density of the 
phase volume (see Eq. (fTT|l '). 

In the classical derivation of Eq. (|20(l we have assumed 
that the spreading of the initially localized energy distri- 
bution during the correlation time Tc is such that the 
corrections to the diffusion coefficient are small in com- 
parison with its zero order value, ^/DtcD' « D. Esti- 
mating the derivative as D' ~ D/E, where is a charac- 
teristic energy of the system, we obtain the condition of 
applicability of the classical equation of energy diffusion, 



E^ 



(21) 



which is much more tolerant than Eq. H12|l . Bridging the 
gap between Eq. I|12|) and Eq. H21|) remains a challenge 
for the quantum theory. 



C. The account of damping 

Now wc turn to the account of damping. For systems in 
which the logarithmic rate of energy damping is constant, 
E = —jE, the energy distribution could be written as 



' {E, t) = 6^*7/. {Ee'^^) , 



(22) 



where 'ip{z) is an arbitrary positive integrable function. 
This functional form can be expressed by the equation in 
partial derivatives 



dQ_ 

dt 



(23) 
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By combining Eqs. ^ and H23() we obtain the equation 
of energy diffusion in the perturbed chaotic system with 
damping: 



= 0. 



fts stationary solution is given by the formulas 



Q{E) 
R{E,Lj) 



7V$ [E) exp[i? {E,io)\ , 

E 

-7 " ^ 







de, 



(24) 

(25) 
(26) 



where N is the normalization constant. The stationary 
distribution could be canonical, Eq. , only in the spe- 
cial case when the diffusion coefficient is exactly propor- 
tional to energy. 



III. THE MODEL 

The equations Eqs. (|25I26() could be checked by di- 
rect numerical solution of the system Eq. ^ for any 
given potential U{r). The calculation of Q{E) is rela- 
tively easy; however, the estimation of the integral in the 
RHS of Eq. H26() is much more demanding, since at ev- 
ery integration step in e the power spectrum Sx(£, (^) has 
to be calculated with sufficient accuracy. Instead of this 
computation we will proceed with analytical approaches 
workable with a special choice of the model. 

We take the PuUen - Edmonds oscillator Q , that de- 
scribes the two-dimensional motion of a particle in the 
quartic potential, as the model for the isolated system. 
The Hamiltonian of this system is 



pI) + 



9 ? 

X y 



(27) 



In the following we use the particle mass m, the frequency 
of small oscillations a-nd the nonlinearity length A as 
unit scales, and write all equations in dimensionless form. 

The properties of chaotic motion of the PuUen - Ed- 
monds model are 

thoroughly studied 0,0, El- With 
the increase of energy the system becomes more chaotic 
both in extensive (that is characterized by the measure of 
the chaotic component /j.^ [E) on the surface of Poincare 
section) and in intensive (that is measured by the magni- 
tude of the Lyapunov exponent <t{E)) aspects. For values 
of energy E > 2.1 the measure fis > 0.5, and chaos dom- 
inates in the phase space; for E > 5 the chaotic motion 
of the system is approximately ergodic p^. 
For the selected model Eqs. © have the form 



X + '-fx + X (l + y'^) = Fsinuit, 
y + jy + y(l + x'^) 0; 



(28) 



the external force is chosen to be directed along the OX 
axis. It may be noted that it is sufficient to couple the 



external force to only one dynamical variable, as opposed 
to the Langevin forces given by Eq. ||2Jl: in chaotic system 
the interaction between vibrational modes will secure the 
redistribution of energy. 

For a particle in a two-dimensional potential the value 
of (f>(i?) is proportional to the area (j){E) of the region 
that is bound by the equipotential line U{x,y) = E, 
namely ^(E) = 2n(j){E). For the PuUen - Edmonds 
model 



/2E 



4>{E)^^ 



2E-x^ 
1 + x^ 



-dx. 



(29) 



This integral can be calculated analytically: 



cj){E) = A\/2E + 1 



K 



2E 



2E 



2E+1 



2E+1 



(30) 

where K(x) and E{x) are complete elliptic integrals of the 
first and second kinds correspondingly. For E ^ 1 this 
function has the asymptotic form 



(E) « V32E [in V32E 1 



(31) 



The approximate expression for the power spectrum 
of coordinate x for the PuUen - Edmonds model was ob- 
tained in Ref. from the assumption of ergodicity of 
motion and in the "frozen frequencies" approximation: 



Sx{E,Lu) 



3(f>{E) 



2E+1 



1 



'1/2 



, (32) 

for 1 < w < V2E + 1 and {E, w) = outside of this 
range. The comparison of this formula with the numeri- 
cally found spectrum is shown in Fig. 1. 



CO 








0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Frequency co 

FIG. 1: The dependence of the power spectrum of the coordi- 
nate Sx of the PuUen - Edmonds model Eq. <\27^ on frequency 
Lo for the energy value E = 16. Theoretical estimate Eq. H32|l 
(line) and numerical calculation (points). 

It can be seen that the main body of spectrum, apart 
from a small interval of frequencies around w = 1, is 
described satisfactorily. 
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IV. THE GIBBS ATTRACTOR 

In this section we study different properties of the mo- 
tion on the Gibbs attractor. For the numerical calcula- 
tion we will mainly use the following set of parameters: 
7 = 2- 10^"^, F = 0.3 and uj = 1.1, that corresponds to 
G = 103. Later it will be referred to as a standard set. 

1. The motion of our model on the Gibbs attractor 
is chaotic: numerical computation gives for its Lyapunov 
exponent ct = 0.65(2). From the Kaplan - Yorke con- 
jecture [20I I21I I it follows that the fractal dimensionality 
of the attractor Dp will differ from the dimensionality 
of the phase space dp = 4 by a quantity of the order 
7/cr = 3-10"'^. This difference is hardly noticeable: from 
the practical point of view the Gibbs attractor densely 
fills the phase space, in resemblance with the canonical 
distribution. 

2. The energy distribution on the Gibbs attractor 
can be determined by substitution of Eq. 132|l into Eq. 
(PS)l : thus we get 



R{E,Lu) = - 



- 1 



£$ (e) 



(2£ + 1 - t^2) 



2^,3/2 



de. (33) 



It is convenient to introduce the temperature parame- 



ter 



e = 



(34) 



that is similar to the temperature of the equilibrium 
created by white noise heat bath and viscous damping, 
Q — D/'j: it scales as the square of the force amplitude 
and as the inverse value of the damping parameter. 

For values of frequency a; > 1, after the substitution 
of the asymptotic form Eq. 1)321) in the numerator and 
disregard of the quantity (cj^ — 1) in the denominator, 
the integration in Eq. (|33|) yields the simple formula: 

R{E,Lu) = -—{lnE + \n32^3). (35) 

In the derivation of Eq. (|24|l we have assumed that 
the energy damping in the absence of the external force 
is governed by the equation E — —^E. For the PuUen - 
Edmonds model this relation must be corrected. The 
exact equation for the energy damping has the form 
E = — 27T, where T is the kinetic energy (properly aver- 
aged), but the virial ratio v{E) — 2T/E in general case 
depends on the energy. From the virial theorem for the 
homogeneous potentials of the power k it follows that 
v{E) — const = 2k/{k + 2). For small energies, £^ <C 1, 
the quartic term in the potential U{x,y) is negligible, 
and v{E) k, 1, but in the domain i? ^ 1, that we are 
interested in, this term becomes important. For purely 
quartic potential we have v{E) = 4/3. Thus we may 
expect that for i? ^ 1 in the PuUen - Edmonds model 
1 < v{E) < 4/3; a naive interpolation leads to value 



v{E) = 7/6. The numerical calculation shows that the 
asymptotic value of v{E) for high energies is very close to 
this value indeed, with the accuracy not worse than 2% 
for £^ > 20. Therefore for very large values of G we can 
improve the expression for R{E, lu) just by multiplying 
it by the factor 7/6; the corrected value of the damping 
parameter will be denoted as 7. 

Comparison of the theoretical form of the probability 
distribution Q{E) (with the virial correction included) 
with its values found numerically is shown in Fig. 2 



0.03 



o 



0.02 




40 60 
Energy E 

FIG. 2: The stationary probability distribution Q{E) for 
the Gibbs attractor of the Pullen - Edmonds model with the 
standard set of parameters (7 = 2- 10^'^, F — 0.3 and uj = 
1.1). Theoretical form given by Eqs. ^E^, and ll^ 

(line) and numerical calculation (points). 

The relaxation time for the diffusion process can be 
estimated as the ratio of the square of the characteristic 
width of the stationary distribution (e.g. the variance 
of the distribution Ve = (E"^) - (Ef) to the average 
value of the diffusion coefficient: r,. = Ve/ {D). For the 
standard set of parameters theoretical estimate based on 
Eqs. Cni), and (EH) gives Ve = 420, (D) = 0.67 
and Tr = 620, whereas the numerical calculation gives 
Ve = 315, (D) = 0.42 and = 740. In our numer- 
ical experiments we have used the time averaging over 
the intervals about 1.7 • 10^ ~ 230rr, that ensures high 
accuracy of the stationary distribution. 

3. In opposition to the canonical distribution, the 
Gibbs attractor is limited in the phase space. From the 
equation for the energy dissipation the average rate of 
energy loss due to damping in the state with the energy 
E is 



P_ = jE. 



(36) 



The maximal rate of increase of the energy due to the in- 
fluence of the external force (for F <C 1) is reached when 
the particle moves along the OX axis with its velocity in 
phase with the external force, and equals to 



F. 



(37) 
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Balancing these quantities we find that the energy of the 
system can not exceed the value 



E. = 



2^ 



(38) 



For the standard set = 8.3 • 10^ = 809. The prob- 
ability to exceed this value with the energy distribution 
given by Eqs. (|25I26|I is about exp(— 750); thus the upper 
limit is not observable. 

On the other hand, the attractor is not limited in the 
energy from below. The basic approximation of nearly 
ergodic motion is not applicable in this range. However, 
the system that at some moment is at rest, in the state 
with E = 0, can be driven by the external force to the 
energies as high as 



Eq 



1 



2(cj2 - ly 



(39) 



(in harmonic approximation) . For the standard set Eq — 
4.72. Thus in our case the efficient communication be- 
tween states with negligibly small energies and those in 
the nearly ergodic domain is possible. 

4. The motion of the system Eq. with the 

standard set of parameters is found to be chaotic up 
to the maximal available times of numerical calculation 
(t = 1.7-10^ for 10 different initial conditions). The time- 
averaged value E of the energy of this motion is -E w 26. 
Changes in parameters {7, F, co} that lead to diminution 
of the average energy of motion on the attractor, such 
as the increase of damping, or lessening of the amplitude 
of the harmonic force, or driving this force out of reso- 
nance with the chaotic unperturbed motion, eventually 
suppress the chaotic motion. For some initial conditions 
the system, after more or less prolonged interval of tran- 
sient chaotic motion with E ^ 1 (we note in passing that 
main results of our theory are applicable to long-living 
transient chaos as well as to the perpetual one), rapidly 
turns to regular motion on a limit cycle or a torus with 
E < 1. This behaviour can have two explanations. 

(A) The Gibbs attractor turns into a semiattracting 
set that supports a metastable chaotic motion. 

(B) The Gibbs attractor endures, but its basin of at- 
traction shrinks, thus ceding larger parts of the phase 
space to basins of regular attractors. 

At present we can not choose between these alterna- 
tives. 

We can establish the border of the domain of the ap- 
parent presence of the Gibbs attractor conventionally, as 
the surface in the parameter space at which for 80% of 
randomly chosen initial conditions the commencement of 
the regular motion takes less than 3 • 10^ units of time. 
The detailed definition of this surface is a laborious task, 
so we have restricted ourselves by variation of each pa- 
rameter in turn, with the other two being fixed at the 
standard values. 

With the increase of damping the domain of presence of 
the Gibbs attractor is limited by the value 7 = 5- 10~^, 



where the transient chaotic motion has E Ki 11. With 
diminishing 7 the chaos persists: even for 7 = the sys- 
tem displays chaotic motion, accompanied by unlimited 
spread of the energy distribution with time. 

With the decrease of the amplitude of the force the 
domain of presence of the Gibbs attractor is limited by 
the value F « 0.11 where the transient chaotic motion 
has E ^ 6. The increase of the force F up to values far 
beyond F ^ 1 does not bring any noticeable changes in 
the character of motion. 

This can be explained as follows. The magnitude of 
the external force F must be compared to the averaged 
(absolute) value of the force in the same direction, cre- 
ated by the static potential, Fu — {\—dU/dx\). On the 
energy surface E it can be calculated by integration: 



/2E 



Ft, 



3<j,{E) 



2E- 



1 + x^ 



2E 



1 



-dx . 



(40) 



For large E this expression has the asymptotics Fu ~ 
E/liiE. The temperature parameter can be written 
as 9 = with k ^ 1, and the average energy of 

chaotic motion in the domain of existence of the Gibbs 
attractor, 9^1, depends on 9 as (E) ~ 9/ln9. 
Thus for the average potential force we have the esti- 
mate (Fu) ^ K^F^ In^^ kF, and for the ratio F/ (Fu) we 
get 



{Fv 



IhkF 
K^F ■ 



(41) 



This quantity is always small and, moreover, decreases 
with growth of F. In physical terms, the harmonic force, 
albeit large in comparison with the scales of the unper- 
turbed model, heats the system to the level of the energy 
content, at which it becomes relatively small in compar- 
ison with the average potential force, thus maintaining 
the energy diffusion picture adequate. 

Finally, on the frequency scale the domain of presence 
of the Gibbs attractor is limited by the band from u « 
0.93 where the transient chaotic motion has i? sa 7, to 
uj w 1.6 where E « 11. 

5. Since the energy of the system in our approach is 
the main dynamical quantity, it is appropriate to study 
its correlation function, 



Be iT)^{EiO)E{T))~{Ey 



(42) 



The first term in the RHS of this equality could be ex- 
panded in the Taylor series in the time shift r and written 
as 



(i?(0)i?(r)) =^if„ — , 



Tl = 



(43) 



where the coefficients Kn are equal to the products of the 
initial energy and of the initial value of energy's n-th time 
derivative, averaged with the stationary distribution. 



E [0)^(0) 



(44) 



7 



From Eq. (24) for the local rate of the energy change we 
have E = -jE + A + D' and 

Ki ^ (eE^ = {--fE"^ + AE + D' E) . (45) 

The higher time derivatives of E can be found by con- 
sequent differentiation of Eq. (24) in time, recurrent 
substitutions and integrations by parts. For example, 

K2 = (eE^, where 

E^{-fE-A- D') (7 - A' - D") + {A" + D'") D. 

(46) 

The normalized correlation function of energy fe^ (r) = 
Be{t)/Be{0), found in the numerical calculation, is 
compared to the theoretical estimate in Fig. 3. 




FIG. 3: The normalized correlation function of energy 6_b (t) 
for the motion on the Gibbs attractor of the PuUen - Edmonds 
model with the standard set of parameters {j — 2 ■ 10"'^, 
F = 0.3 and u — 1.1). The time shift r is measured in 
periods of harmonic field. Theoretical approximation by first 
three terms of the expansion Eq. 1431 (line) and numerical 
calculation (points). 

The rate of decay of the energy correlations is deter- 
mined by a competition of the absorption, that domi- 
nates in the lower energy range, and of the loss of energy 
through damping, that prevails in the upper range; they 

are balanced by the stationarity equation (e\ — 0. The 



dependence of this rate on the parameters of the model 
is rather complicated. 



V. CONCLUSION 

The study of the system of the type Eq. ||SJ) may be 
important, since for many experimentally relevant Hamil- 
tonian models one can indicate some mechanism of relax- 
ation, that could be approximated by the viscous damp- 
ing terms at least qualitatively. 

The obvious candidates for the further studies of the 
Gibbs attractors are chaotic billiards. Our equations per- 
mit us to get some simple estimates for this case. Since 
a billiard per se has only two dimensional parameters, 
the particle mass m and some characteristic size a, the 
scale of time depends on the initial conditions and is pro- 
portional to E~^/'^. Thus the power spectrum of the co- 
ordinate can be written (with m and a as units) in the 
scaling form 

where g{z) is some positive integrablc function. From 
discontinuities of velocity for the law of motion x{t) it 
follows that the high frequency asymptotics of the power 
spectrum has the form Sx{E,uj) oc w""*. Combining this 
formula with Eqs. (EHJ, and (gZ)), we obtain the 
approximate form of the exponent in the stationary dis- 
tribution: 

R{E,uj)^-^, (48) 

where Q = cF'^j^^uj^'^ and c is a numerical constant that 
depends on the specific form of the billiard. This estimate 
is valid for E <^ uj^ for billiards of any dimensionality. 
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